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We have developed an ensemble density functional theory which includes spin degrees of freedom 
for nonuniform quantum Hall systems. We have applied this theory using a local-spin-density ap- 
OS ' proximation to study the edge reconstruction of parabolically confined quantum dots. For a Zeeman 

OS \ splitting below a certain critical value, the edge of completely polarized maximum density droplet 

reconstructs into a spin-unpolarized structure. For larger Zeeman splittings, the edge remains po- 
<—j ' larized and develops an exchange hole. 



I. INTRODUCTION 

The recent development of nanofabrication technology has made it possible to manufacture semiconductor systems 
with reduced dimensionality and very high electron mobility. These technological advances-have led to the discovery 
J) ' °f such fascinating phenomena as the fractional and integer quantum Hall effects (QHK-jEJ. These occur in a two- 
dimensional electron gas in a magnetic field B = Bz perpendicular to the electron systemcl. A quantum treatment of 
the motion of an infinite, homogeneous system shows that the kinetic energy takes discrete values (n+ l/2)tko c , where 
n is the Landau level index (n = 0, 1, 2, . . .) and uj c — eB/m*c is the cyclotron frequency. Each Landau level contains 
Cj ng = B/$q states per unit area, or one state for each magnetic flux quantum <j> = hc/e, giving rise to a macroscopic 
I ' Landau level degeneracy. The ratio of the electron areal density n(r) to tie defines the filling factor z/(r) = n{v)/riB- 
The filling factor also can be expressed as v — 2irl B n, where Ib = \jficjeB is the magnetic length. 

The fractional quantum Hall effect can occur when electron-electron interactions dominate disorder. At certain 
filling factors of the form v = p/q, with p and q relative primes and q odd, electron-electron interactions cause 
the condensation of the electrons into highly correlated states. These states are incompressible with an energy gap 
separating the ground state from the bulk excited states. However, in a finite system, there must exist gapless 
excitations localized near the edgescl. Thus, the low-energy physics of finite systems is dominated by the gapless 
edge modes. Therefore, in is necessary to be able to accurately model edges of FQHE systems in order to explain 
experiments. In a finite FQHE system with the potential confining the electrons varying slowly compared to lg, 
the electronic structure 9-trthe edges may form a series of alternating compressible and incompressible regions with 
t-H ■ a step-like density profilcO'tl In addition to standard transport measurements, there are now a variety of probes to 
directly study edge structures in inhomogeneous systems. Examples are capacitance spectroscopy of the quantum 
Hall edgeaJ, time-resolved measurements of edge magnetoplasmonsQ, and surface acoustic waves techniques which are 
capable of resolving very small spatial inhomogeneities in the electron densityH. Addition spectroscopy has also been 
used to study quantum dots with sizes of the order of 100 nm and with 10 to 100 electronsti 

For an explanation of experimental studies it is highly desirable to have a computational approach which accurately 
treats systems with of the order of 1 — 10 3 electrons, and which can include effects such as accurate confinements, 
i spin degrees nf freedom, and finite layer .thickness. Exact numerical diagonalizations are limited to very small systems 
(N < 1^,)c3'EJ. Semiclassical methodsO'El do not accurately treat electron-electron interactions, and effective field 
theoriesL3 cannot give accurate quantitative information about many system properties. A method that can deal 
with larger number of electrons is the composite fermion theory in the Hartree approximationE 3 ! However, in this 
approach, the singular Chern-Simons gauge field is replaced by its smooth spatial average, and the composite fermion 
mass has to be put into the calculations by hand. Furthermore, interpretation of the results is sometimes difficult 
and ambiguous. On the other hand, density functional theory (DFT) is known as a general quantitative method to 
include exchange-correlation effects in inhomogeneous systems without any fitting parameters. In this paper we show 
that it can be used to give highly accurate results for quantum Hall systems. Preliminary results for spin polarized 
systems were reported earlierE-l'ILa 

The DFT was originally formulated J&y Hohenberg and Kohn as a practical method for a description of the ground 
state properties of many-body systemsta . The foundation of DFT is the Hohenberg-Kohn theorem, which states that 
the ground state density uniquely determines the Hamiltonian of a system. Furthermore, a variational principle states 
that ±he ground state density minimizes the energy of the system. We will use the constrained search formulation of 
LevyO for the Hohenberg-Kohn theorem and its associated variational principle. In this elegant approach the ground 
state energy E can be written as a functional of density 



o 
o 



> 

OS 



X 



1 



E[n] = F[n] + J dr n(r)F ext (r). (1) 

Here 

F[n] = inf (*|f + Vee|*>, (2) 

*— >n 

with T, 14e, and V ex t kinetic energy, electron-electron interactions, and external potential, respectively. The infimum 
is taken over all many-body states \t that yield a fixed density n(r). F[n] so defined is then a universal functional 
of the density n(r). For a given external potential V cx t, the true ground state density is the function n(r) which 
minimizes E[n] in Eq. (|l|). ,— ■ 

The origins of the DFT are to be found in the statistical method developed by Thomas and FermitS. They first 
realized the advantage of describing an inhomogeneous systems by using the density 

u(r) = J dr 2 . . . J dr N | *(r, r 2 , . . . , rjy) | 2 , (3) 

an observable, rather than the unobservable complex wave function '5 of Nd variables in d dimensions. Thomas- Fermi 
theory is a way to find an approximate n(r\__The theory is valid in the semi-classical limit and has a successful history 
of applications to many different problems&J. But since Thomas-Fermi method neglects exchange-correlation effects 
and makes an approximation for the kinetic energy functional, it has some serious deficiences as well. For example, 
Thomas-Fermi theory cannot predict ferromagnetism. DFT remedies these problems by explicitly (and formally 
exactly) incorporating exchange-correlation effects as well as interaction parts of the kinetic energy functional into an 
exchange and correlation energy functional E xc [n(r)], and by developing a useful computational scheme for including 
exchange-correlation effectscl This is done by introducing an auxiliary non-interacting system with a ground-state 
density n s (r), and by asserting that there exists an effective potential V s (r) for this system such that n s (r) — n(r), with 
n(r) the ground-state density of the real, interacting system. A system with this property is called w-representablc. 
The density is then obtained from a simple Slater-determinant of the so-called Kohn-ShamfKS) orbitals i/j a (r), 
n s (r) = ^2a=i IV'a ( r ) 1 2 j where - 0a( r ) are obtained by self-consistently solving the KS equationsEJ 

/i eff Va (r) = [T + V s (r)] ^ a (r) = s a ijj a (r) . (4) 

The self-consistency is achieved in practice by iteratively obtaining the eigenstates and occupying the N eigenstates 
with the lowest eigenvalues e a . The effective potential V s (r) can be derived from the Hohenberg-Kohn theorem and 
its associated variational principle. First, the functional E[n] [Eq. (|l|)] can be decomposed as 



E[n] = T [n] + / drn(r) 



V ext (r) + ^Vu(r) 



(5) 



Here To is the kinetic energy of a non-interacting system with density n, Vh is the classical (Hartree) Coulomb 
potential 

with eo the static dielectric constant, and E xc may be viewed as a definition of the exchange-correlation energy. The 
non-interacting kinetic energy T is treated exactly in this approach, which removes many of the deficiencies of the 
Thomas-Fermi model. The variational principle applied to Eq. (|^) yields 

S jn^ = ^h + ^ext(r) + Yk(r) + V KC (r) = M , (7) 
on(r) on(r) 

where [i is the Lagrange multiplier associated with the requirement of constant particle number, and the exchange- 
correlation potential V xc is formally defined as the functional derivative 

. , SE-xdn] , „ 

V m (r) = -=^U. 8 
on{r) 

Comparison of Eq. (R) with the corresponding relationship for a non-interacting system, 
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on[r) on(r) 

gives an expression for the effective potential 

V s (r) = V cxt (r) + V a (r) + V xc (r). (fO) 

In Eq. (|J) all terms but the exchange-correlatiop-energy E xc can be evaluated exactly. In practical calculations, the 
local density approximation (LDA) is often usedcJ. In this approximation, one writes 

£xc° A = / drn(r)e xc (n(r)), (11) 

where e xc (n) is the exchange-correlation energy per particle of an infinite, homogeneous system of density n. The 
exchange-correlation potential [Eq. (||)] is then obtained as 



,ldA/„a d[ne xc (n)} 



dn 



(12) 

,(r) 



Above we ignored the electron spin - n(r) was the total electron density, and the spin degree of freedom was ne- 
glected. The Hohenberg-Kohn theorem formally ensures that every property, including the spin density or polarization, 
can be obtained from the ground-state density. However, practical LDA calculations of systems with spontaneously 
broken symmetries, such as spin rotation symmetry, typically are much improved if the order parameter of the broken 
symmetry, e.g., spin density or polarization, is explicitly included by construction. In particular, the broken symmetry 
may not otherwise be obtained accurately from the LDA. In GaAs samples, where most of QHE experiments have 
been done, the spin degree of freedom is important, and may lead to inhomogeneous spin densities. This is because 
the effective Zceman energy g* fisB is quite small compared to the cyclotron energy fiu) c , g* (J-bB jhio c k, 0.02, due 
to the small effective mass m* = 0.068m e and a reduced Lande g-factor \g*\ — 0.44. This means that in a uniform 
noninteracting system, two highly degenerate Zeeman levels with the same Landau level index n are almost degenerate 
in energy. Since the dielectric constant for GaAs eo ~ 13, the cyclotron energy Tiw c and Coulomb energy e 2 /(eo?s) 
are of the same order of magnitude for magnetic fields of the strength of a few Tesla ( ^/("o^) ~ 0Ay/B[T]). Using 
the Coulomb energy as the unit of energy, the dimensionless parameter characterizing Zeeman coupling can then be 
defined as 



(13) 



- = 9*VbB 
9 ~ e 2 /Ms)' 

The small value of g (typically about 0.02) makes the existence of partly polarized states energetically possible^ 
even at v < 1. Therefore, it is necessary to include the spin degree of freedom in any quantitative theoretical 
approach to QHE systems. A reasonable first step,is to generalize the DFT to include the spin polarizationEj'Ei. An 
exact treatment of the spins, in general, requiresE3 the replacement of the charge density n(r) by the single-particle 
density matrix p aa >(r) — (0\ip£ (r)Vv(r)|0). Here, ifi£ (r) and Vv(r) are the usual field operators corresponding to the 
annihilation and creation of an electron with spin a at r, and |0) is the ground state of the system. With a constant 
magnetic field applied in the z-direction, the z component of the total spin angular momentum is a constant of the 
motion and it is convenient to assume that the magnetization density only has a z component. Under this assumption, 
the single-particle density matrix can be taken toJae diagonal, p aa i (r) — p aa i (r)6 aa i . In this case the constrained 
search procedure of Eqs. ([l]) and (|2|) is modified tcc3 

£[n T ,nj] = F[n h n L ] + J drn(r)U cxt (r), (14) 

where F[n|,nJ = inf^^^lT 1 + Ez + Kd'I') and Ez is the Zeeman energy with ^ yielding fixed densities n a (r) . 
The local spin-density approximation (LSDA) is then given by 

E^ DA [n hni ] = J drn(r)e xc [n r (r), ni (r)}, (15) 

where e xc [nj,rtj is the exchange-correlation energy per particle in a homogeneous system with up- and down- spin 
densities and nj,, respectively. In spite of the LSDA being justified only in the limit of small spatial variations of the 
electron density, this approximation has been surprisingly successful in describing the properties of inhomogeneous 
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atomic, molecular and solid-state syst This scheme correctly predicted, for example, ferromagnetism in Fe, 

Co and Ni among the transition metalst£l. Moreover, the self-interaction-corrected LSDA was successfully applied to 
some strongly correlated systems. such as the transition-metal oxides and a Hubbard model representing a Cu02 layer 
in the cuprate superconductorsEZI. 

In the original Kohn-Sham formulation [Eqs. (|l|) and (0)], the ground state was assumed to be nondegenerate and the 
ground state density n(r) assumed to be pure-state w-representable. This means that n(r) can be be expressed in terms 
of a single Slater determinant of KS orbitals i/j a {r) obeying an effective single-particle Schrodinger equation [Eq. 
However, there are systems which-are known-aot to be u-representable in this sense. One class of such systems was 
considered independently by LevyE3 and LiebEj. Consider a system with q independent iV-particle degenerate ground 
states l^i), . . . , \ ty q )- Then construct the density matrix 

D = J2d^ l )(^ l \ (16) 

i=i 

with di = dl > 0, X)?=i <k = 1- This yields the density 

9 

n(r) = Tr{I)n(r)} = ^ d;n;(r), (17) 

i=l 

where Uj(r) = (5 , i|n(r)|5'i) (spin degrees of freedom are neglected for simplicity). As was proven by Levy and Lieb, if 
q > 2 the density n(r) cannot be represented by a single ground state in DFT, i.e., it cannot be obtained by a single 
Slater determinant of the N lowest energy KS orbitals. However, there exists a generalization of Hohenberg-Kohn 
theorem which provides a one-to-one correspondence between a ground state density n(r) and the Hamiltonian even 
for system with a ground state density which can be of the form of Eq. (0). By extending the functional F[n] in 
Eq. (|) to 

F E [n] = inf Tr{£>(f + V ee )}, (18) 

£>— vn 

with the infimum taken over all D yielding a fixed density n(r), there is then a generalized variational principle which 
states that FeI^] is minimized by the ground state density, which can now be represented by an ensemble of wave 
functions, even if it cannot be represented by a single Slater determinant. This generalization is called ensemble 
density functional theory. As we shall see below, fractional QHE systems are not w-representable, so an ensemble 
DFT has to be used. 

In Sec. H we review our ensemble DFT scheme. The essential features of our ensemble DFT approach will be 



illustrated in Sec. [II by applying it to spin-polarized quantum Hall dot. Sec. [V then describes the spin-unpolarized 
edge reconstruction of v — 1 quantum Hall dot. The phase diagram for the spin textured edge reconstruction of the 
maximum density droplet will be presented there. Finally, conclusions are given in Sec. 0. 



II. PRACTICAL ALGORITHM FOR ENSEMBLE DENSITY FUNCTIONAL THEORY 

In practical ensemble DFT calculations one introduces as in the KS scheme an auxiliary non-interacting system 
which provides the basis for the density matrix and has a ground state density identical to_the interacting system 
at hand. By using the variational principle, one then arrives at a set of equations analogously to the KS equations 
Eq. (|J). However, the density for N electrons is now given by 

n(r)=^/ a |^ Q (r)| 2 , ]T f a = N, (19) 

a a 

with the occupation numbers f a in the interval < f a < 1. One obtains fractional occupancies f a only when the 
corresponding KS eigenvalues e a are degenerate and equal to the Fermy energy ep. (If s a < sf, then f a = 1.) Let us 
show briefly why applying DFT to the FQHE inevitably requires ensemble DFT. We consider an infinite, homogeneous 
fractional QHE system at a filling factor of v = 1/3, and assume that we have the exact exchange-correlation potential 
Vxc for this system. We then construct a set of determinantal iV-particle wavefunctions {^i} made up from the KS 
orbitals ipa obtained by solving the KS equation using the exact exchange-correlation potential. Label each member 
of this set by the set of numbers {8i a } identifying the KS orbitals used to construct ^ [0*™ = (0)1 for (un)occupied 
orbitals, and J2a @ia = N]. Because the KS orbitals are orthogonal, a state ^ has densityEil 
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rii(r) =^0 M |^ Q (r)| 2 . (20) 

a 

Construct three states i = 1,2, 3, with occupancies {6> ict } = {100100 . . .}, {010010 . . .}, {001001 . . .}. The occupied 
KS orbitals are all degenerate lowest Landau level states (since the system is homogeneous), so the three ^ are 
degenerate. Then construct a density matrix of the form of Eq. ( |l6| ) with d\ = d-2 = da = 1/3. The resulting density 
n(r) (Eq. |l7| ) is constant: n(r) — ^ ^ , corresponding to a uniform system with v = 1/3. Now we can appeal to 

the Levy-Lieb theoremUS: since the v = 1/3 ground state density can be constructed as in Eqs. @ and (0) with 
q > 2 degenerate states "J^ this density is therefore not pure state i>-representable, and ensemble DFT must be used. 

In inhomogeneous FQHE systems, not all KS orbitals are degenerate, but some are. By a simple extension of the 
argument in the paragraph above, such systems too are not pure state w-representable, and ensemble DFT must be 
used. For inhomogeneous systems one finds M orbitals with e a < Sf and D degenerate orbitals with e a = £p. One 
constructs determinantal wavefunctions in which all M low-energy orbitals are occupied; the ^ differ by which 
N — M of the D degenerate orbitals are occupied. Using Eq. (|l7|) and the density of determinantal wavefunctions 
given in Eq. (^0|), the total density for the ensemble represented by D can be calculated as 

i 

n(r) = ^^a a |i(r)| 2 . (21) 

a i=l 

Comparing the result with Eq.(|l^), one can see how the fractional occupational numbers f a of the degenerate KS 
orbitals follow from the weights di in the expansion of density matrix D: 

i 

1=1 

As we mentioned in Sec. [l| for the density defined by Eq. (|2l]), a generalization of Hohenberg-Kohn theorem exists and 
an extended variational principle [Eq. (|l8|)] can be used. However, a procedure to compute the fractional occupancies 
f a has not existedo, and one mai o!: . a dvance in our work is that we have found a simple way to generate the 
occupancies, at least for the FQHFxirE.3. Applying ensemble DFT to the FQHE, we have found that fractionally- 
occupied KS orbitals are indeed degenerate at the Fermi energy, consistent with our demonstration above that the 
FQHE is in general not pure-state i>-representable. We will review this scheme here. 

In our algorithm, we start with a set of input occupancies and single-particle orbitals and iterate the system N eq 
times using the KS scheme. The number N cq is chosen large enough (about 20-40 in practical calculations) that the 
density is close to the final density after N eq iterations. If the density of the system could be represented by a single 
Slater determinant of the KS orbitals, we would now essentially be done. However, in an "ensemble w-representable" 
(but not pure state w-representable) system there are many degenerate or near-degenerate KS orbitals at the Fermi 
energy, and small fluctuations in the density between iterations cause the KS scheme to occupy a different subset of 
these orbitals each iteration. This corresponds to constructing different Slater determinants each iteration. In 
other words, when the KS orbitals are degenerate at the Fermi energy there is an ambiguity in how to occupy these 
degenerate orbitals. Nevertheless, by associating the number of the iteration j with the particular Fock state 
represented by the occupation numbers 8j a , one can accumulate the running average occupancies after each iteration 

where N- lt is the number of iterations. This is clearly of the form of Eq. j2^ ) with the weights di given by the relative 
frequency of the corresponding Slater determinant during N it — N cq iterations. To achieve a self-consistent solution, 
the system of KS equations is then iterated using average occupancies f a ,N it instead of 8j a , with the density after 7Vj t 
iterations calculated as 

nN it (r) = ^2f a , N Mr)\ 2 - (24) 

a 

This density is used to calculate the new effective potential V s (r) according Eq. (]To|) , and by solving the resulting 
KS equations a new set of occupancies 8j a is obtained, which can be again used to update averages f a ,N it ■ When a 
self-consistent solution is obtained this highly nonlinear map should give the same input and output densities after 
one iteration. To prove formally that this scheme with the running average occupancies converges self-consistently to 
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the physical density is by no means a trivial problem. However, we have numerically verified for FQHE systems that a 
finite-temperature version oi-pur scheme converges to a thermal ensemble at finite temperatures down to temperatures 
of the order of 10~ 4 hu! c /kB^- We have also performed Monte Carlo simulations about the ensemble obtained by our 
scheme. In these simulations, we used a Metropolis algorithm to randomly change the occupation numbers about 
our converged solution, keeping the chemical potential fixed. The free energy of the new set of occupation numbers 
was calculated self-consistently. The results were that to within numerical accuracy our scheme gives the lowest free 
energy. We have also checked our ensemble DFT against small system numerical diagonalizations with very good 
results (see Sec. ED). 

This algorithm has made ensemble DFT a practical calculation tool, and it may be possible to apply it to strongly 
correlated systems other than the QHE, for example inhomogeneous Mott insulators, provided there are accurate 
approximations for the ground state energies of homogeneous systems available. 



III. SPIN-POLARIZED QUANTUM DOT IN A FQHE REGIME 

In order to illustrate the essential features of our ensemble DFT scheme, we first neglect the spin degree of freedom 
and consider a spin-polarized quantum dot in the FQHE regime. Typical quantum dots contain about 10-100 electrons, 
so the quantum dot can be used as a model system in order to demonstrate the usefulness of the DFT in the study of 
large inhomogeneous electron systems. Moreover, quantum dots are believed to have highly correlated ground states 
in strong magnetic fieldsa. Hence, these systems can be used ptOj-show how well these strong correlation effects are 
represented by LDA compared to exact diagonalization studieal3'tll. There have been some attempts to use density 
functional theory to model such systems. Ferconi and Vignaleo performed current-spin DFT studies of quantum dots 
in the integer QHE regime with a small number of electrons. In their calculations, the energy gaps due to correlation 
effects were not included, and the spin-dependence of the exchange-correlation energy was taken to be that of an 
electron gas in the B — > limit. Their results for a spin-polarized three-electron system are in good agreement with 
results from exact diagonalizations. However, their approach cannot be extended to include fractionally occupied 
states or the complicated spin dependence of the exchange-correlation energy in the strong magnetic field region. 

Electron-electron interactions in the ground state of quantumjdots have been studied experimentally, for example by 
measuring the tunnelingrconductance through a Coulomb island&E3. The dots used experimentally can very accurately 
be modeled as parabolicQ, i.e., the Coulomb islands are confined by a parabolic potential V ex t — m*Q 2 r 2 /2, where m* 
is the effective mass of an electron, characterizes the strength of confining potential and r is the distance from the 
center of the dot. We chose this potential for our model system. 

Due to the circular symmetry we can label the KS orbitals ip a (r) by Landau level index n > 0, and by angular 
momentum label m > — n (a = {m, n}) and expand ip a (r) in the eigenstates \mn) of the single-particle Hamiltonian 
To = (-ihV + f A) 2 /2m* as 

in(r) = e""V mn (r) = ^C mnn >(r\mn'). (25) 

n' 

We can then solve self-consistently the KS equations h e sip mn — e mn (p mn for the radial parts <p mn of KS orbitals for 
each value of the angular momentum m separately. The orthonormal basis (r\mn) can be written in the cylindrical 
gauge A = ^B(x, —y, 0) in terms of the associated Laguerre polynomials^ L™ 



(r\mn) 



L™(^- T )e^-r A /^. (26) 



These are centered on circles of radii r rnn w \j2(m + ti)Ib with Gaussian fall-offs for r <C r mn and r 3> r mn . We 
write the density n(r) in dimensionless form as a nonuniform filling factor v(r) = 2i:l 2 B n(r), 

V = 2nl B S ' / W m|V , mra| 2 = ^ ' fmnC m nn' C mrLn " (mn 1 77171 ). (27) 



mnn n 



This is used to calculate the effective potential V s (r) in the KS equations. In the lowest Landau level these expressions 
simplify (for example, the filling factor becomes just v(r) = e~ x J2 m f m0 l^T w ith x = r 2 /4Z^). However, we have 
used the four lowest Landau levels (n = 0, . . . , 3) and kept the general expressions, since as we mentioned in Sec. Q, 
the cyclotron energy %lo c and Coulomb energy e 2 /(eo^s) are of the same order of magnitude for typical experimental 
situation, so there can be appreciable Landau- level mixing. The solution of the KS equations in the basis \mn) reduces 
to iteratively obtaining the eigenstates with the lowest energies of the following block-diagonal Hamiltonian matrix 
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h e ff,mn'n = hu c {n + -) + (run \V s \mn) , (28) 

where the effective potential V s (r) = V cx t(r) + Vn(r) + V xc (r) is defined by Eq.(|Io|). We chose a confining potential 
of strength Ti£l = 1.6 meV corresponding to that of McEuen et alB. The Hartree potential for a circularly symmetric 
quantum dot is given by 

Va(r) = -4- d <t> I dr> ; (29) 
eo's Jo J V^ 2 + r' 2 — 2rr' cos(/> + <5z 2 

where is a finite thickness of the 2D electron liquid. The interaction of the 2D electrons with a gate can also be 
easily incorporated into the Hartree potential. Although the typical value of Sz for a real experimental situation is 
about 10 nm, in our simulations we took 6z — > in order to compare our results with other theoretical works. 

We have used the LDA to obtain the exchange-correlation potential V xc . In this approximation, we first need the 
exchange-correlation energy per particle in a homogeneous system e xc which we chose as 

e xc M = ^rV) + 4M- (30) 
The first term is a smooth interpolation formula of Lcvesque et 

eS m M= r*r(-^-Wr)-l] 
Jo \reohJ 

~ -0. 782133^ (1 - 0.211i/ 0,74 + O.OUu 1 - 7 ) (e 2 /e Z B ) (31) 

for the ground-state energy obtained by evaluating the pair correlation functions g v {r) at certain fillings v < \ for 
about 256 particles using very accurate Monte Carlo methods. The second term in Eq. (p0|) , e^ c (v), contains the cusps 
in the ground state energy which cause the FQHE. The discontinuity in the slope of e xc {v) near certain "magic" filling 
factors v* = p/q is. related to the chemical potential gap A^i — q(\A p \ + |A^|). Here A Pj /j are the quasiparticle (hole) 
creation energiestS at v = v* . In our calculations, we restrict ourselves to include only the cusps at v — 1/3, 2/5, 3/5 
and v = 2/3, which are the strongest fractions. (See Appendix A for a detailed description of our expression for 
e£..) Substituting Eq. (|3^) into Eq. ( |l2|) we can find the exchange-correlation potential as a function of filling factor 
V xc (v) for a spin-polarized system. Thispotential is depicted in Fig. (H) as V xc< i(v, £ = 1) (the generalized spin DFT 



approach will be discussed later in Sec. |V|) 



The discontinuities in V xc {v) in the LDA give rise to a numerical instability. The reason is that an arbitrarily 
small fluctuation in charge density close to an FQHE fraction gives rise to a finite change in energy that leads to 
serious convergence problems. To overcome this, we made the compressibility of the system finite, but very small, 
corresponding to a finite, but very large, curvature instead of a point-like cusp in e xc at the FQHE fractions. This 
was accomplished by allowing the discontinuity in chemical potential to occur over an interval in filling factor S of 
magnitude 10 -3 . This corresponds to a sound velocity of about 10 6 m/s in the electron gas, which is three orders of 
magnitude larger than the Fermi velocity of a 2D electron gas at densities typical for the FQHE. In general, the finite 
compressibility does not lead to any spurious physical effects so long as the energy of density fluctuations on a size 
of the order of the systems size is larger than any other relevant energy in the problem. The only noticeable effect 
is that incompressible plateaus, at which the density would be perfectly constant were the compressibility zero, have 
density fluctuations on a scale of S. 

We have used our ensemble DFT-LDA scheme to study the edge reconstruction of a 40-electron quantum dot as, 
a function of magnetic field strength. For a certain range of magnetic field (typically about 2 - 3 T), it is knowncEl 
that such a dot forms a so-called maximum density droplet. The maximum density droplet is spin-polarized with 
a filling factor which is unity in the bulk and falls off rapidly to zero at tq w \J2NIb- It is the most compact 
droplet (minimum angular momentum) that can be formed of spin-polarized electrons in the lowest Landau level. 
Increasing the magnetic field increases the importance of electron repulsion compared to confinement. At some 
higher magnetic field B, the edge of maximum density droplet is reconstructed, forming an exchange hole, because 
it becomes energetically favorable, to spread out the electron density while still taking advantage of short-ranged 
attractive exchange interactionocJ. 

For a, -slowly varying confining potential, alternating compressible and incompressible strips in density may be 
formedu'u between v = 1 and v = regions, with the density of the incompressible strips fixed at the density of 
an FQHE fraction, v = p/q. The widths of the incompressible strips are determined by the energy gaps at the 
fractions p/q, while the widths of the compressible strips are determined by electrostatics. These compressible and 
incompressible strips at the edge of an FQHE system were studied in an extended Thomas-Fermi approximation at 
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low but finite temperatures by Ferconi, Geller, and Vignalec^l for infinite Hall bars. Their results for the widths of the 
incompressible and compressible strips were in good agreement with the predictions by the semiclassical theoriesutl. 
In the extended Thomas-Fermi approximation, the kinetic energy was treated as a local functional, as in the standard 
Thomas-Fermi approximation, while the exchange-correlation energy was included in a LDA. This extended Thomas- 
Fermi approximation is presumably valid in the limit of very slowly varying confining potential and large numbers of 
electrons. In contrast, our ensemble DFT approach (Sec. ||) treats the kinetic energy exactly, and does not have any 
limitations on the number of the electrons. Furthermore our calculations can be done at very low or zero temperature. 

Figure [l] depicts some of our results for a 40-electron droplet using our ensemble density functional approach. In 
these calculations, we chose a parabolic confinement of strength hil — 1.6 meV, and eo = 12.4, appropriate for GaAs, 
and T = 0. For values of the magnetic less than about 2.5 T, the droplet forms a maximum density droplet. AtBw 2.8 
T, an exchange-hole forms. This value of magnetic field compares very well with values obtained from Hartree-Fock 
calculations. For example, using eo = 12-4 and hfl — 1.6 meV, we obtain from the work of MacDonald, Yang, and 
JohnsonH a value of 2.52 T, while the value from the work of de Chamon and WenO is 2.84 T. As the magnetic 
field is increased further, the effective confinement softens, and the droplet undergoes several reconstructions in which 
compressible and incompressible strips are formed. Figure ^ shows one example of incompressible strips at v = 2/3 
and v = 3/5 for a magnetic field of 4.1 T. In our calculations, the incompressible strips form at filling factors where V xc 
has a discontinuity. Traversing across an incompressible strip, the effective single-particle potential (V eK t + Vh + V xc ) 
is constant - as V ex t + Vh varies, V xc varies across its discontinuity exactly to screen the change in V^xt + Vh- Thus, 
the width of the incompressible strips is determined by the distance over which V eyLt + Vh varies by .-an amount equal 
to a discontinuity in V xc , in agreement with the argument of Chklovskii, Shklovskii, and Glazmara. The widths of 
the compressible strips are determined by electrostatics plus the smooth parts of the exchange-correlation potential. 

In order to compare our results with exact diagonalization studies of quantum Hall dots by Yang, MacDonald and 
Johnsont3 we have calculated the expectation value of the total angular momentum (M) = ^2 mn mf mn as a function 
of the magnetic field strength B for N = 6 spin-polarized electrons in the lowest Landau level and a confining potential 
of hfl = 2.0 meV. The results are shown in Fig. ||for two different versions of the exchange-correlation energy e xc . The 
diamonds (o) were generated using the Levesque-Weiss-MacDonald exchange-correlation energyE3, while the plusses 
(+) were generated using an exchange-correlation energy due to Fano and Ortolanie3. This latter is constructed not 
only to give the right behavior for v < 1/2, but both the value of e xc at v = 1/2 and particle-hole symmetry of the 
lowest Landau level are included to give a good interpolation formula on the entire interval < v < 1 for electrons in 
the lowest Landau level. Both exchange-correlation energies give clear plateaus or plateau-like structures in angular 
momentum vs. magnetic field. However, the Levesque-Weiss-MacDonald is a rather poor approximation near v = 1 
(a region for which it was not constructed), and furthermore overestimates the magnitude of the exchange-correlation 
potential at about v = 1/2. As a consequence, the initial maximum density droplet instability is smeared out and 
starts at a too low value of magnetic field, and the formation of a 1/3 droplet (as is evidenced by studies of the density 
profile) occurs at a too high value of magnetic field. Also, the values of the angular momentum at the plateau-like 
regions tend to be too low. For example, the formation of the 1/3-droplet occurs at M w 40, while the exact value is 
M = 45. In contrast, the results obtained using Fano-Ortolani exchange-correlation energy tend to be very accurate. 
For example, the maximum density droplet instability occurs at B ss 2.8 T in our calculations, compared to B — 2.75 
T in the numerical diagonalizations, and the 1/3 droplet formation occurs at B ss 5.3 T in our calculations, compared 
to B — 5.29 T in the numerical diagonalizations (see Fig. ||). In addition, the plateau- like regions are more developed 
and flatter in angular momentum. Still, though, the ensemble DFT tends to underestimate the angular momentum 
at the plateaus. We want to emphasize here that we have not used any adjusting parameters in our calculations. 
Furthermore, the ensemble DFT is not constructed so as to give only integer angular momentum. Finally, only the 
energy gaps at v = 1/3 and v = 2/5, along with their particle-hole conjugates at v = 2/3 and v = 3/5 were included, 
while the numerical diagonalizations used the full Coulomb interaction. Therefore, perfect agreement between our 
ensemble DFT results and the numerical diagonalizations cannot be expected. We also did these calculations at a 
finite temperatute of 100 mK (ksT / 'e 2 / (cqIb) ~ 1 x 10 -3 ), which improves the convergence of these small particle 
systems. In view of all this, the agreement between our ensemble DFT results and the numerical diagonalizations 
of Yang, MacDonald and Johnsont2l must be considered remarkable. We are presently working on extending the 
Fano-Ortolani interpolation to include v > 1 and several Landau levels. 



IV. SPIN TEXTURED EDGE RECONSTRUCTION OF THE MAXIMUM DENSITY DROPLET. 



The edge reconstruction of v = 1 maximum density droplet outlined in Sec. Ill suggests that v — 1 is much simpler 
to study than the fractional quantum Hall filling factors due to the absence of a cusp in e xc (v) at v = 1. Therefore, 
one might expect that v = 1 is a Fermi liquid in the sense that the elementary excitations are well described by single- 
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particle excitations and only renormalized by the interactions. However, recent experiments on high mobility GaAs 
quantum wellsCj have provided evidence for the existence of topological charge-spin textures, so-called 'skyrmions', 
near v = 1. These are non-trivial many-body excitations due to electron-electron interactions first predicted to be the 
low-energy excitations near v — 1 by Sondhi et alc3, with the energies about half of those of single-particle spin-flip 
excitations. The fact that skyrmions are the low-energy excitations near v — f (and also possibly near v = 1/3) 
raises the possibility of spin-textured edge reconstruction of the maximum density droplet. Therefore, inclusion of the 
spin degree of freedom may be essential in the study of inhomogencous systems. Indeed, Hartrcc-Fock and effective 
field theoretical calculations have shown that for a soft confining potential, the edge of an infinite Hall bar at v = 1 
becomes unstable to spin-textuxed reconstruction for weak Zeeman coupling, while stronger Zceman coupling yields 
a spin-polarized reconstructiono. 

Motivated by these ideas, we have generalized our ensemble DFT approach to include spin degrees of freedom 
within the LSDA [Eq. Jl5[)]. We have applied this generalization to parabolic maximum density droplets to study 
the edge reconstruction as a function of confinement strength (magnetic field) and Zceman coupling. We find that 
as the Zeeman coupling is decreased, the edge becomes unstable-to spin-textured reconstruction at a strength of the 
Zeeman coupling consistent with that found by Karlhede et a/.c3 This provides further evidence for the importance 
of spin degrees of freedom in inhomogeneous systems and demonstrates the usefulness of our LSDA ensemble density 
functional approach. 

For a parabolic dot, the variational principle applied to the KS functional [Eq. (|l4|)] yields two sets of KS equations 
for spin up and spin down electrons 

(r), (32) 

where 



V B , a (r,B) = ag*^B + V cxt (r) + V H (r) + V xc , a (r, B) (33) 
is an effective potential for the auxiliary non-interacting system. In the LSDA the exchange-correlation potentials are 



d 

V XC!<T (r,B) = ^— (ne xc [n T , n il B]) 



(34) 

n a =n a (r) 



The parametric dependence on the magnetic field B can be incorporated by using spin filling factors v a — 2nl 2 B n a as 
variables instead of spin densities n a . To make connection with the spin-polarized case we first transform the spin 
filling factors v a to total filling factor v and spin polarization £ : 



(35) 



v = v\ + v i 

The exchange-correlation potentials (Eq.([}4|)) then become 

d d 

Kc,T = tt" (^e xc ) + (1 - O-^^xc , 
ov ot; 

Vic,i = ^-(^ xc )-(l+0^exc , (36) 

where the exchange-correlation energy per particle in a homogeneous system with a filling factor v and polarization 
£, i.e., e xc = Exc^, £); nas to be approximated. 

We already have given the expression for the exchange-correlation energy per particle for the spin-polarized electron 
gas e xc (z/, £ = 1) [Eq. (p0|)]. The question is then how to obtain a reasonable interpolation formula for e xc between 
spin-polarized (£ = 1) and spin-unpolarized (£ = 0) 2D electron liquids for a fixed v in a strong magnetic field. We 
have constructed a reasonable first approximation, as we now explain. 

In what follows x(c) as a subscript denotes exchange (correlation) respectively. We decompose the exchange- 
correlation energy _B XC into exchange E x and correlation E c energies. Since the exchange interaction only acts between 
parallel spins, we have 

ExWi^i] = ^h.^tl + \ E y\ v l^i\- (37) 

Moreover, it follows from dimensional analysis that the exchange energy must scale as density (filling factor) to the 
3/2 power in a 2D electron gas. Following Oliver and PerdewcH, we can then write: 







/ d A r [^ /2 (r)+"?» 

We also have from Eqs. ( |35| ) i/f = 5^(1 — 0; v i = 1^(1 + 0- Equation ( |38| ) can then be rewritten as 

£ x ~ / dW/ 2 [(l+£) 3/2 + (l-£) 3/2 )" . 



3/2/ 



(38) 



(39) 



Since in the local density approximation for the exchange energy E K = J d 2 rve x (u, £) we are then led to the form 



e*(i>, = exfo e = 1) + (e*(f , C = 0) - e x (i/, £ = 
= e x (i/,C = l)+Ac x (i/,0, 



where the function 



f(0 = 



(l + 3/2 + (l-Q 3/2 -272 
2 - 2^2 



(40) 



(41) 



is an interpolation function between the two extreme cases £ = and £ = 1 with /(0) = 1 and /(l) = 0. Although 
the analogous simple closed form for the correlation energy e c (v,£) is not available, it can be always be written as 
£ c(^, = £ c(i / , £ = 1) + Ae c (z/, £). So, as a first approximation we will use the form of Eq. ( |40| ) for the smooth part 
of the correlation energy e c , too (leaving the cusps aside for the moment), with the same interpolation function /(£), 
as was suggested first by von Barth and HedinE3. Denoting the smooth part of the exchange-correlation energy per 
particle by e xc we can then write 



>,0 



i) 
i) 



[4 

5e x 



>,e=i)]/(0 



(42) 



So far, we have constructed a function e xc (f, £) which gives a smooth interpolation for the exchange-correlation 
energy for any value of v and £. What is left is to add the cusps to this function. We already have a good approximation 
for these at £ = 1. We now need to extend this approximation to arbitrary values of £. Very little is known about 
the cusps, i.e., the energy gaps, for arbitrary polarizations. It is known that there is a gap for un-polarized systems 
at fillings v = 2/5, v = 3/5, and v = 2/3. The gap, and thus the cusps, occur at very special 'magic' configurations 
at which the system can take advantage of a particularly low correlation energy. Therefore, it seems plausible that 
for a given value of v, say v = 2/5, there cannot be an energy gap for any value of £ between and 1. In order to 
incorporate this assumption into a usable approximation, we interpolate our cusp energy constructed for polarized 
systems, e xc (i / ), to arbitrary polarizations by multiplying it by a function g(£) which is unity at £ = and £ = 1 with 
zero derivative at these points, vanishes away from these values of polarization, and is symmetric about £ = 1/2. All 
together, then, we have 



iWM 



Specifically, we chose 



9(0= [4£ 2 -l] 2 [27-£ 2 (40 - 16£ 2 )] /27, 



(43) 



(44) 



which is the only polynomial in £ satisfying the above constraints. Near v = 1 (where there is no cusp in the total 
exchange-correlation energy), the sign of the function <5e xc (^) will then determine the spin polarized (ferromagnetic) or 
spin unpolarized (paramagnetic) ground state of the infinite electron liquid (neglecting the Zecman splitting). Indeed, 
substitution of Eq. ( [i"3| ) into Eq. (^) gives 



XC,| 



xc, J. 



25e xc (v)f'(0 + 2eg(v)g'(0. 



(45) 



The last term in this expression may be ignored near v = 1. We would thus expect a spin-polarized ground state if 
<5exc(^) > because /'(£) < for all < £ < 1 so in this case the inequality V xc ,t < V*c,i holds. Otherwise we would 
expect a spin unpolarized state. Since a von Barth-Hedin type approximation has been aunlied before only to the 
3D electron gas in zero magnetic field, it is natural to ask how faithful this approximation iscj to data obtained from 
the numerical diagonalization studiesE3c3 outlined in Table Q. First, reversed spins are in fact possible in the ground 
state of quantum Hall liquids. It is known from numerical studiesEJ that the electron gas is completely polarized only 
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at filling factors v = 1, ^, Secondly, there are some partially polarized ground states that cannot be realized in the 
approximation used here because the function /(£) is monotonic. If the partially polarized states were realized by this 
function it would have a minimum at some fractional £. However, the simple model Eq. (42) allows us to capture the 
essential physics of the spin unpolarized edge reconstruction of the quantum dot as it will be shown on the example 
below. 

In order to obtain the function 8e^ c {v) in Eq. ( |43"|) wc start by calculating the energy differences between spin 
polarized and unpolarized states for fractions listed in Table [j]. The value for the ground state energy of a v = 1 
unpolarized system is not available, but a reasonable approximation is to take e xc (^ = 1, £ = Q) = e xc (i> = 1/2, £ = 1) = 
—0.469 (e 2 /eo^s), implying that the spin-up and spin-down components are uncorrelated. The ground state energy of 
a v = 1 polarized system e xc (i/ = 1,£ = 1) = —0.6265 (e 2 /e /s). Therefore we have <5e xc (V = 1) = 0.1575 (e 2 /c Ib)- 
at quantity is close to the energy required for a exchange-enhanced single spin flip, y/ V / 8(e 2 / e$l g) , at this 
To complete the numerical parameterization of exchange-correlation functional, we then perform a spline 
fit to obtain the function <5e xc (^). We have plotted the exchange-correlation potentials V KCy(7 as a function of a filling 
factor v at £ = 1 and £ = in Fig. ||. We see that at £ = 1 the difference between exchange-correlation potentials 
for spin up and spin down electrons AV XC [Eq. ([45])] is changing sign from negative to positive while the filling factor 
v is decreasing from v = 1 to v — 2/3. Ignoring the Zeeman splitting and the cusps, the ground state of an infinite 
electron liquid would change from spin polarized to spin unpolarized. To estimate the possibility of having a spin 
unpolarized state above filling 2/3 with the inclusion of the Zeeman splitting, we have to compare the dimensionless 
Zeeman energy g [Eq. @] with the difference AV XC [Eq. @] at this filling AV X cO = 2 /3,£ = 1) » 0.05(e 2 /e Z B ). 
This value is larger then the Zeeman splitting for GaAs g « 0.02. Therefore, the ground state of a GaAs based 
homogeneous system is a spin unpolarized state at and just above filling factor 2/3. In an inhomogeneous system, 
in addition to exchange-correlation potential and Zeeman energy, there are also the Hartree interaction of the 2D 
electrons and the external potential which confines them. Hence, even in von Barth and Hedin type approximation, 
it is possible to have not only polarized and unpolarized states, but also a partially polarized state in inhomogeneous 
system such a quantum dot. 

We have investigated the spin density of a quantum dot using our LSDA ensemble DFT. The results, choosing the 



same parameters as in Sec. Ill, are shown in Fig. ^[ We find that the electron liquid is indeed partially polarized 
at the quantum dot edge provided the Zeeman energy is small enough (Fig. For Zeeman energies below a 
certain critical value, g < g c , the maximum density droplet is reconstructed with the increasing of the magnetic 
field by forming an unpolarized (in general partially polarized) state. However, as g increases above the critical 
value, g > g c , the polarized reconstruction described in Sec. Ill takes place. We can plot the phase diagram of 
the edge reconstruction of the quantum dot in the (3,7) plane (Fig. where the dimensionless parameter 7 = 
m*Vl 2 l B / '(e 2 / 'cqIb) characterizes the 'softness' of the edge. The phase boundaries separate the maximum density 
droplet, spin polarized and unpolarized instabilities. The value of g c ~ 0.055 separating the spin-polarized and spin- 
structured instabilities is in good agreement with- the value gJ MJ ~ 0.03 found from numerical diagonalization of 
parabolic dots byr-¥ang, MacDonald and JohnsonEB. A phase diagram analogous to ours was obtained in the work by 
A. Karlhede et alE3. Their value of the critical Zeeman splitting g^ KLS = 0.169 obtained by Hartree-Fock calculations 
is about 3 times larger then g c w 0.055 from our phase diagram. We speculate that the difference is due to the the 
fact that correlations were ignored in their Hartree-Fock calculation, and to the different geometry they used - an 
infinite Hall bar. 



V. SUMMARY AND CONCLUSIONS 



We have developed a spin ensemble density functional approach to strongly correlated systems and used it to 
study inhomogeneous quantum Hall systems in the integer and fractional Hall regimes. For spin-polarized systems, 
our approach gives results in excellent agreement with numerical diagonalizations, Hartree-Fock and semi-classical 
calculations. Note that while all of these latter approaches have limited regions of applicability, such as small systems, 
systems near v = 1, or the semi-classical limit, we have here demonstrated that our ensemble density functional 
approach spans all these regions, which makes it a useful approach to general inhomogeneous quantum Hall systems. 

We have generalized the ensemble DFT to include spin degrees of freedom within a simple local spin density 
approximation, and applied this generalization to a quantum dot. Our results show that for small, but physical, 
Zeeman energies, g < g c , the maximum density droplet is unstable with respect to spin-textured edge reconstructions 
as the magnetic field increased. At larger Zeeman splittings, g > g c , the maximum density droplet is unstable 
with respect to spin-polarized edge reconstructions. Our value of g c is in good agreement with that obtained from 
numerical diagonalization studiestlJ. Hartree-Fock calculations for an infinite Hall bar by A. Karlhede et a/.c3 give a 
phase diagram qualitatively analogous to ours. However, Hartree-Fock calculations are limited to v ps 1, while our 
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ensemble DFT is in principle applicable to general fractional quantum Hall systems, e.g. droplets at v = h. The 
accuracy of our approach depends on obtaining good estimates of the exchange-correlation energy as a function of 
both electron density and spin polarization for homogeneous fractional quantum Hall systems. Work is currently in 
progress to obtain such estimates. Finally, the spin ensemble DFT used here cannot be used to study spin-charge 
textures (skyrmions), in which the spin polarization rotates smoothly in space. Work is currently under way to 
generalize our spin DFT to include such charge-spin textures. 
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discussion about Ref. |28[ and J.M. Kinaret for help with the LSDA. O.H. would like to thank Chalmers Institute of 
Technology, where part of the work was done. This work was supported by the NSF through grants DMR93-01433 
and DMR96-32141. 



APPENDIX 



We will construct the cusp-part of the exchange-correlation energy, e xc (i/) for a spin-polarized system by first 
considering v < 1/2, and then use electron-hole symmetry to obtain the form for 1/2 < v < 1. Finally, for v > 1 we 
assume e xc (z^) = e xc (l — v). 

For spin-polarized systems in the lowest Landau level, we write e xc (z/) e xc (;/) + t^ c (v), where e xc (i/) 1S given by 
a smooth interpolation, such as the Levesque-Weiss-MacDonald formulaEZI (although this one does not obey strict 
particle-hole symmetry in the lowest Landau level), or the Fano-Ortolani formulaEa Particle-hole symmetry yields 
for the total exchange-correlation energy 

v Mi/) - e xc (l)] = [1 - u] Ml -v)- e xc (l)] , (46) 

from which we obtain 

i/O) = i/'O*) + (1 - 2^)e xc (l), (47) 

with v* = 1 — v. This means that 

i/O) = (48) 

We define 

g(u) ee i/4(i/). (49) 

Since the discontinuities in the chemical potential at fractional QHE fillings p/q is a relation for d[ve xc (v)] /dv, it is 
easier to work with g(v) than with e xc (i/). Then particle-hole symmetry implies that 

dj, = _d_g^) 

av dv* 

At fractional QHE fillings p/q, we must have 

^ 0(*OU(p/,)+ " ^ ff(")U(p/ g) - = + (51) 

where /i + and /j_ arc the quasi-particlc and quasi- hole creation energies (defined to be positive), respectively. We 
construct g{v) for < v < 1/2 to be piece- wise smooth, with <?(^) = for p/q a fractional QHE filling, and a 
discontinuity in the derivative given by Eq. (^l|) . 

We only included the cusps at v = 1/3, 2/5, their particle-hole conjugates, and the corresponding values at fillings 
inreased by unity. For v < 1/3, we make the Ansatz 

g{v) = avq^ (p/q) (y - p/q) |>-PAd - 3o ] , (52) 

with p/q = 1/3 and 

a = W^o) (53) 
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and g = exp(-p/q). 

For 2/5 < v < 1/2 we take 



0(1/) = -M2/5) [l-e^" 2 /*)] , (54) 



a 



with ao = 80. 

Next, for 1/3 < v < 2/5 we used a cubic interpolation 

g{u) = a{v - 1/3) (i/ - 2/5)(i/ - i*). (55) 
Fixing the slope of g(z/) at (l/3) + and (2/5) - then yields 

„ . 3/*+(l/3) 



(1/3 -2/5) (1/3-^)' 

and 



(56) 



_ 5M2/5)/3 + 6 M+ (l/3)/5 ^ 
M+ (l/3) + 5^(2/5) ■ 1 } 

Finally, we smooth out the resulting discontinuities in V xc over an interval 26 about the fractional QHE fillings 
v = p/q. To do this, we interpolate linearly g{v) between its values at v = p/q±5, so that dg(y)j dv = A+B(is—p/q+5), 
where A = g(p/q — 5) and B = [g(p/q + S) — g{p/q — 6)] /(2<5). Simple integration then yields g{v) = Av + \Bv 2 — 
Bv(£ — S) + C, where C is an integration constant given by 

C = g( P /q -S)- A(p/q -S) + \B{pjq - S) 2 . (58) 
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Potential Energy 
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7 


-0.3870 




-0 3868 
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Unpolarized 
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-0.4410 
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-0.0057 


Unpolarized 
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-0.3975 
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-0.3970 
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Partially polarized 




-0.4219 
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-0.4241 


-0.0022 


Partially polarized 


4- 


-0.4528 
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Partially polarized 




-0.5232 
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-0.0099 


Unpolarized 
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-0.5010 


-0.5044 




-0.5096 


-0.5074 


0.0064 


Partially polarized 



TABLE I. Potential energy (per particle) for various values of spin polarization and 5t xc {v) (Eq.( ^2|)) for the four-electron 
(rows 1-6, Ref.[20]) and six-electron (rows 7 and 8, Ref. [46]) systems. The Zeeman energy is not included in the energy values. 
The unit of energy is (e 2 /coZb). 
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FIG. 1. Edge reconstruction of a spin-polarized quantum dot as the magnetic field strength is increased. Plotted here is the 
local filling factor v{r) for a parabolic quantum dot with TiQ. = 1.6 meV and 40 electrons. For magnetic field strengths B < 2.5 
T the dot forms a maximum density droplet, and for B ~ 2.8 T, an exchange hole is formed. For stronger magnetic fields, 
incompressible regions form, separated by compressible strips. 

FIG. 2. Expectation value of the total angular momentum (M) = ^ mf mn as a function of the magnetic field strength B 
indicated by solid line for a spin-polarized six-electron droplet in a parabolic confinement using the Levesque-Weiss-MacDonald 
(Ref. 37) (<>) and Fano-Ortolani (Ref. 42) (+) exchange-corrleation energies. The solid shows the exact diagonalization studies 
result from Ref. 10. 

FIG. 3. Local filling factor as a function of radial coordinate for a six-particle system in a parabolic external potential with 
HQ = 2.0 meV. Here, the Fano-Ortolani exchange-correlation energy (Ref. 42) was used. The transition to a 1/3-droplet occurs 
between B — 5.3 T and B = 5.4 T. The bumb in electron charge at the edge of the system is characteristic of systems with a 
not too soft confining potential. 

FIG. 4. The exchange-correlation potentials V xc ,<t as a function of a filling factor v at £ = 1 and £ = in units of e 2 /(cpls)- 
The solid line indicates Vx C ,t an d short-dashed line corresponds to V xc ,i at £ = 1. According to Eqs. (^) and (|4l|), the 
exchange-correlation potentials V^ c ,a coincide at £ = (since /'(0) = 0) and are shown by the long-dashed line. The increase 
in V^ c as functions of v at a FQHE filling factors occurs over a range of a filling factor of 0.002. 



FIG. 5. Spin structured instability at the edge of a quantum dot for a Zeeman splitting g = 0.014, magnetic field B — 3.05 T, 
and N = 38. The external potential is characterized by TiQ. = 1.6 meV (so that the dimensionless strength of the confinement 
is 7 = 0.063). (a) The solid line depicts the total local filling factor v(r) as a function of radial coordinate r, and the dashed 
line depicts the polarization, (b) The occupancies of the KS states (/i m o,o(r) are plotted against orbital center coordinate 
r m = (2m) 1//2 ?B with "+" for majority (|) spin occupancies, and "o" for minority (j) spin occupancies. At the instability of 
the maximum density droplet for this value of g, there is a minority-spin population at the edge of the system, (c) Eigenvalues 
of the two lowest Landau level KS orbitals, with + depicting eigenvalues of the majority spin orbitals, and o depicting the 
eigenvalues of the minority spin orbitals. The chemical potential is indicated by the solid line. At the edge the filling factor 
takes fractional values, and the KS eigenvalues are here degenerate and equal to the Fermi energy sf, in agreement with the 
general theory of Sec. II. 

FIG. 6. Phase diagram of the edge reconstruction of a parabolic quantum dot in the (g, 7) plane for N = 38 electrons. Here, 
the confining potential has a strength given by hQ, = 1.6 meV. For 7 > 0.065, the system forms a maximum density droplet for 
all values of the Zeeman coupling g. For values of the Zeeman coupling g larger than a critical value g c , the maximum density 
droplet undergoes an initial reconstruction to a spin-polarized exchange hole as the confinement strength 7 is decreased, while 
for g < g c the maximum density droplet has a spin-structured instability with decreasing 7. In these calculations, g c ~ 0.055. 
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